The temperature dependence of the electron mobility in molecular crystals

We show that a hopping model may explain the experimental temperature dependence of the electron mobility in molecular crystals. The flat, higher temperature portion is obtained if anharmonic effects are included in a simple manner. The low temperature part is obtained if coherent effects are added to the diffusional ones.     

Inelastic scattering with Chebyshev polynomials and preconditioned conjugate gradient minimization

We describe and test an implementation, using a basis set of Chebyshev polynomials, of a variational method for solving scattering problems in quantum mechanics. This minimum error method (MEM) determines the wave function Psi by minimizing the least-squares error in the function (H Psi - E Psi), where E is the desired scattering energy. We compare the MEM to an alternative, the Kohn variational principle (KVP), by solving the Secrest-Johnson model of two-dimensional inelastic scattering, which has been studied previously using the KVP and for which other numerical solutions are available. We use a conjugate gradient (CG) method to minimize the error, and by preconditioning the CG search, we are able to greatly reduce the number of iterations necessary; the method is thus faster and more stable than a matrix inversion, as is required in the KVP. Also, we avoid errors due to scattering off of the boundaries, which presents substantial problems for other methods, by matching the wave function in the interaction region to the correct asymptotic states at the specified energy; the use of Chebyshev polynomials allows this boundary condition to be implemented accurately. The use of Chebyshev polynomials allows for a rapid and accurate evaluation of the kinetic energy. This basis set is as efficient as plane waves but does not impose an artificial periodicity on the system. There are problems in surface science and molecular electronics which cannot be solved if periodicity is imposed, and the Chebyshev basis set is a good alternative in such situations.     

Density functional study of the charge on Aun clusters (n=1-7) supported on a partially reduced rutile TiO2(110): are all clusters negatively charged?

It is widely believed that small gold clusters supported on an oxide surface and adsorbed at the site of an oxygen vacancy are negatively charged. It has been suggested that this negative charge helps a gold cluster adsorb oxygen and weakens the O-O bond to make oxidation reactions more efficient. Given the fact that an oxygen vacancy is electron rich and that Au is a very electronegative element, the assumption that the Au cluster will take electron density from the vacancy is plausible. However, the density functional calculations presented here show that the situation is more complicated. The authors have used the Bader method to examine the charge redistribution when a Aun cluster (n=1-7) binds next to or at an oxygen vacancy on rutile TiO2(110). For the lowest energy isomers they find that Au1 and Au3 are negatively charged, Au5 and Au7 are positively charged, and Au2, Au4, and Au6 exchange practically no charge. The behavior of the Aun isomers having the second-lowest energy is also unexpected. Au2, Au3, Au5, and Au7 are negatively charged upon adsorption and very little charge is transferred when Au4 and Au6 are adsorbed. These observations can be explained in terms of the overlap between the frontier molecular orbitals of the gold cluster and the eigenstates of the support. Aun with even n becomes negatively charged when the lowest unoccupied molecular orbital has a lobe pointing in the direction of the oxygen vacancy or towards a fivefold coordinated Ti (5c-Ti) located in the surface layer; otherwise it stays neutral. Aun with odd n becomes negatively charged when the singly occupied molecular orbital has a lobe pointing in the direction of a 5c-Ti located at the vacancy site or in the surface layer, otherwise it donates electron density into the conduction band of rutile TiO2(110) becoming positively charged.     

Does phenomenological kinetics provide an adequate description of heterogeneous catalytic reactions?

Phenomenological kinetics (PK) is widely used in the study of the reaction rates in heterogeneous catalysis, and it is an important aid in reactor design. PK makes simplifying assumptions: It neglects the role of fluctuations, assumes that there is no correlation between the locations of the reactants on the surface, and considers the reacting mixture to be an ideal solution. In this article we test to what extent these assumptions damage the theory. In practice the PK rate equations are used by adjusting the rate constants to fit the results of the experiments. However, there are numerous examples where a mechanism fitted the data and was shown later to be erroneous or where two mutually exclusive mechanisms fitted well the same set of data. Because of this, we compare the PK equations to "computer experiments" that use kinetic Monte Carlo (kMC) simulations. Unlike in real experiments, in kMC the structure of the surface, the reaction mechanism, and the rate constants are known. Therefore, any discrepancy between PK and kMC must be attributed to an intrinsic failure of PK. We find that the results obtained by solving the PK equations and those obtained from kMC, while using the same rate constants and the same reactions, do not agree. Moreover, when we vary the rate constants in the PK model to fit the turnover frequencies produced by kMC, we find that the fit is not adequate and that the rate constants that give the best fit are very different from the rate constants used in kMC. The discrepancy between PK and kMC for the model of CO oxidation used here is surprising since the kMC model contains no lateral interactions that would make the coverage of the reactants spatially inhomogeneous. Nevertheless, such inhomogeneities are created by the interplay between the rate of adsorption, of desorption, and of vacancy creation by the chemical reactions.     

Rotatable anisotropy of epitaxial Fe<Subscript>1−x</Subscript>Ga<Subscript>x</Subscript> thin films

We show by a combined magnetic force microscopy and synchrotron radiation spectroscopy study that stripe-like patterned magnetic domains are present in Fe_1?x Ga _x thin films. These stripes, whose origin is attributed to an out-of-plane magnetic component, can be rotated by an external magnetic field.     

The boundary layer in the hub area of wind turbine blades at high wind speeds

At low tip speed ratios, the strongest effect at the inboard half-span, when a vortex like flow occurs, is produced. These 3-D rotational effects contribute to a stall delay as well as higher lift and angles of attack at the inboard sections, as compared to 2-D flow conditions. A greater understanding of the mechanisms and properties of the shear layer close to the upper side surface of blade is sought through the development of an analytic theory and then is completed by CFD computations. The results show that the secondary flow has a slip effect on the blade surface and its circumferential velocity is less than that of blade revolution, particularly important for TSR <3.0. While at the high tip speed ratios the flow around the blade is much as the 2-D one, then at the low tip speed ratios, TSR < 3.0 the flow behaves rather like the impinging flow on a rotating disc. Therefore, the 3-D rotational effects dominate the flow at high wind speeds and a special approach of such flows would be to assume. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The change of the vibrational frequencies of a chemisorbed diatomic caused by the electrostatic interactions with the metal surface

Results of density functional calculations are used to discuss the change of the vibrational frequency of an adsorbed molecule caused by the electrostatic interaction between a chemisorbed molecule — simulated by two effective charges — and a metallic surface.     

A one-dimensional microscopic model for the rate of thermal desorption of an atom. The role of multiphonon processes

A one-dimensional microscopic model is used to calculate the desorption rate and to investigate the role of multiphonon processes in desorption. If the vibrational frequency of the chemisorptive bond is larger than the Debye frequency, multiphonon processes of all orders must be included as in our one-dimensional microscopic multiphonon model.